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Motivated by the recent discovery of superconductivity in Ca- and Yb- intercalated graphite (CaC6 
and YbCe) and from the ongoing debate on the nature and role of the interlayer state in this class of 
compounds, in this work we critically study the electron-phonon properties of a simple model based 
on primitive graphite. We show that this model captures an essential feature of the electron-phonon 
properties of the Graphite Intercalation Compounds (GICs), namely, the existence of a strong 
dormant electron-phonon interaction between interlayer and tt* electrons, for which we provide a 
simple geometrical explanation in terms of NMTO Wannier-like functions. Our findings correct the 
oversimplified view that nearly-free-electron states cannot interact with the surrounding lattice, and 
explain the empirical correlation between the filling of the interlayer band and the occurrence of 
superconductivity in Graphite-Intercalation Compounds. 

PACS numbers: 74.25.Jb, 63.20.Kr, 74.70.Ad, 63.20.Dj 



Introduction 

Superconductivity in Ca- and Yb-intercalated graphite 
(CaCg and YbCg), with T c of respectively 11.5 and 
6.5 K, was discovered^ five years later than in the 
pseudo-graphite MgB2, where T C =40K. 3 Superconduc- 
tivity in MgB 2 is understood to be caused by the strong 
interaction between a few zone-centered holes in the 
B sp 2 cr-bonding band and a few zone-centered bond- 
stretching phonons, which themselves get softened by the 
interaction^^ 7 - the tt bands are also at the Fermi level, 
but they couple to phonons only weakly. What makes 
MgB2 special is that the electrostatic and covalent attrac- 
tion of the B p z tt electrons by the intercalated Mg ++ ions 
is so strong that the center of the tt bands is lowered below 
the top of the a band, whereby ZioZe-doping takes place 
into this bonding band. An equally simple and elegant 
understanding of the origin of superconductivity in the 
strongly electron-doped Ca.Ce and YbCg is still missing. 
Csanyi et al& have empirically related the appearance of 
superconductivity in Graphite-Intercalation Compounds 
(GICs) to the filling of a free-electron-like interlayer (IL) 
band, 9,10 which is present, albeit empty, in pure graphite; 
and proposed that plasmon-mediated superconductiv- 
ity could take place in these compounds. Recent exper- 
iment o 11 ' 12 and ab-initio calculation o 13 : 14 ' 15 : 16 agree, in- 
stead, that the conventional e-ph coupling suffices to ex- 
plain the superconductivity in CaCg and YbCg, although 
some puzzling issues remain^ 

It seems therefore that in graphite, if the Fermi sur- 
face has only tt sheets, the e-ph coupling is weak and 
T c < IK, whereas, if also a or IL sheets are present, the 
e-ph coupling increases enough as to yield much higher 
T c 's. Mazin et alJ£ and Calandra et al^ argued that 
the main role in the e-ph coupling is played by electronic 
and vibrational states associated with the presence of 



the intercalant ion, and, generally speaking, this is true: 
without the intercalant, graphite is not a superconduc- 
tor. However, a close inspection of the e-ph spectral 
functions a 2 F(iv) of CaCe of Refs. [T3fl5l shows that a 
considerable fraction (60%) of the total e-ph coupling 
comes from modes which involve carbon atoms and not 
the intercalant ion. This observation, together with the 
empirical correlation between the filling of the IL band 
and superconductivity, suggests that, in analogy with the 
hole-doped graphite compound MgB 2 , CaC6 and YbCe 
might be understood as "electron-doped graphite" , where 
a different source of e-ph interaction is switched on when 
electrons are driven into the IL band. 

In this paper we will show, via density-functional- 
theory calculations^ that indeed a strong interband e-ph 
interaction between tt* and IL electrons, due to out-of- 
plane phonons, survives in graphite even when the in- 
tercalant atom is replaced by a uniform positive back- 
ground of charge (J ellium- Intercalated- Graphite, JIG). S 
Such a finding contrasts with the rough idea that the IL 
states, being free-electron like, cannot interact with the 
surrounding lattice: this may become true when the car- 
bon layers are very far apart from each other, but not 
at the spacings observed in real GICs, where, as our re- 
sults clearly show, the IL states are strongly affected both 
by the presence and by the dynamics of the carbon lay- 
ers which sandwich them. Similar confined, free-electron 
states are found in several layered material o 19 i 20 ' 21 i 22 and 
even in nanotubes^ 3 - where our results suggest that sim- 
ilar e-ph coupling interactions may take place. 

The paper is organized as follows: in the first section 
we discuss the electronic structure of JIG with the aid of 
Wannier functions. In the second and third sections we 
calculate the vibrational spectrum and e-ph interaction of 
JIG. We then consider these as a function of doping and 
interlayer spacing in the fourth section, comparing the re- 
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suits of our model with the available experimental data 
on GICs. Finally, in Appendix A, we provide a minimal 
3x3 tight-binding model for JIG, while Appendix B con- 
tains the computational details for the different methods 
used: planewaves and pseudopotentials, 24 TB-LMTO 25 
and NMTO.— 



I. ELECTRONS 
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FIG. 1: (Color Online) Band structure of primitive graphite. 
The fatness (and color) of the bands is proportional to the 
partial contribution of each orbital channel: green (a), red 
(ir) and blue (interlayer = IL).— The lattice constants used 
in this plot are those of CaC6 (see text and note 28). The 
Fermi level of CaC6 would in this simple picture be at about 
4.5 eV, see also Fig. [9] 

A good starting point to understand the more com- 
plicated band structure of actual GICs, such as CaC6 
and YbC6, is represented by primitive graphite. In this 
structure (space group 191) the carbon hexagons sit on 
top of each other (aaa stacking), i.e., they are not lat- 
erally shifted like in real graphite; the primitive cell thus 
contains two carbon atoms (C2). 

The color and fatness^ of the electronic bands of C2, 
shown in Fig. [T] indicate the relative importance of the 
partial contribution of each orbital channel to the given 
band: green (a), red (ir), and blue (IL). 

With aaa stacking and negligible dispersion perpen- 
dicular to the sheets, the bonding ir and antibond- 
ing 7T* bands touch at a (Dirac) cone in energy-versus- 
momentum space; with an occupancy of 4 electrons per 
C, the Fermi surface consists of merely the K-points in 
the Brillouin zone; simple graphite is an ideal semimetal. 
In the simplest rigid-band picture, hole-doping shifts the 
Fermi level down, towards and eventually inside the a 
bands, like in MgB2. Electron-doping shifts the Fermi 
level up, and eventually brings it into the IL band (CaCg, 
YbC 6 , and other GICs). 

The main difference between the rigid-band picture 
and actual GICs is that the intercalation of positive ions 



between the graphene layers shifts the IL bands down in 
energy with respect to the n bands, without substantially 
modifying their shape in the relevant energy windowj&i£ 
this effect is reproduced quite well even if one replaces 
the intercalant ions by a uniform background of positive 
charge (JIG), provided that the interlayer spacing is fixed 
at the corresponding value of the actual GIC4 

While the nature of the a and ir bands has been 
extensively discussed in literature, less information is 
available about the IL band in GICs. IL states have 
mainly been described as three-dimensional, nearly-free 
electrons in the empty interlayer region, whence their 
name; and indeed, as far as the motion parallel to the 
graphite planes is concerned (the xy directions), the IL 
electrons, to a very good approximation, freely propagate 
in the interstitial volume, with plane-wave-like wavefunc- 
tion exp(ik x x + ik y y) and energy dispersion \(k\ + fc 2 ). 

In the z direction, however, their motion is not free, be- 
ing periodically damped, with period c, by the graphene 
layers, across which the IL electrons are only partially 
transmitted. Along z, these layers act like the periodic 
potential barriers of the Kronig-Penney model, 29 where 
the free-electron parabola \k\ is split into a set of sub- 
bands separated by gaps, with sub-bandwidths and gaps 
dictated by the strength of barriers and by their spa- 
tial separation c. In this model all gaps close, and the 
full 3D free-electron dispersion is recovered, for vanish- 
ing barrier strength; at fixed c, as the barrier strength 
increases, the transmission across potential barriers de- 
creases, until, for infinite strength, electrons are confined 
between pairs of barriers, yielding a set of completely 
flat sub-bands; in our case this translates into IL bands 
£ m = ^(fc 2 +fc 2 ) + i(m7r/c) 2 ,?7i = 1,2, ... with no disper- 
sion along k z . 

A careful examination of different GICs (and of C2 for 
several values of c/a) confirms that their IL bands stay 
somewhere between these two limiting cases, with the 
graphene layers approximately behaving, along z, as an 
array of barriers with fixed finite strength and variable 
period c, which depends on the intercalant atom. These 
bands have a substantial dispersion along k z -so there is 
significant transmission across the graphene layers, or, in 
other words, the IL electrons are not completely confined; 
but their slope and width strongly depend on c -so the 
partial reflection against the graphene layers is significant 
too. 

In the range of parameters analyzed in the present 
paper, which correspond to existing donor-intercalated 
graphites, the lowest m — 1 Kronig-Penney band is well 
separated in energy from the others -another evidence 
that, in the z direction, the nearly- free-electron theory 
does not hold for the IL states- and it is the only one 
which can be occupied; to this band we refer, in the fol- 
lowing, as the IL band (blue in Fig.[T]). 

In Fig. [2] we show the shape of its Wannier function^. 
Because the IL band overlaps other bands -the 7r* and 
a* bands, with which it, however, hybridizes very little,— 
and because the function shown in Fig.[2]has not yet been 
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FIG. 2: (Color Online) Top and side view of the Wannier- 
like^ function V ) s( r )i see a l so Appendix A, for the IL band 
(blue in Fig(TJ . The basis set spanning the IL band con- 
sists of the replicas of this function placed at all interstitial 
sites (grey and blue spheres). The figure shows the contour 
|^ 8 (r)| = 0.07x (a.u.)' 3 , with / iplofr = I. The sign of tp s (r) 
is indicated by red and green. 



symmetrically orthogonalized to the lattice-translated 
functions, it ought to be called Wannier-like as in Refs. 
[H and The basis set spanning the IL band, interpo- 
lated smoothly across any hybridization gap, consists of 
this Wannier-like function and its replicas displaced by 
all lattice translations. We have chosen the IL Wannier- 
like function to be centered on an interstitial site, where 
it is seen to have the full point-group symmetry and a 
simple, s-like shape. The bond-centered Wannier-like 
function for the tt* band, 31 interpolated smoothly across 
hybridization gaps, is the antibonding linear combination 
of the two Cp z orbitals at the vertices of the bond, plus 
a halo of further antibonding C p z orbitals. The basis set 
spanning the tt* band consists of this Wannier-like func- 
tion and its replicas displaced by all lattice translations. 
That is, every second C-C bond carries such a function. 
Since the IL Wannier-like function has everywhere the 
same sign, while the sign of the tt* Wannier-like func- 
tion alternates within the scale of the IL function (there 
is one IL function for every two C p z orbitals), the two 
functions, one displaced against the other by an arbitrary 
lattice translation, are very nearly orthogonal. This is the 
reason why the IL and 7r* bands hardly hybridize, as seen 
in Fig. [TJ On the contrary, the Wannier-like function for 
the 7T band^i is the bonding linear combination of the 
two Cp z orbitals at the vertices of the bond, and that 
is not orthogonal to the IL function. For that reason, 
the tt band may hybridize with the IL band, as will be 
discussed in Appendix A. This hybridization shows up if, 
by doping or changes in the c/a ratio, the IL band moves 
into the energy region of the tt band. 

Based on the band and orbital picture of primitive 
graphite (C2) discussed so far, in the following sections 
we will study the electron-phonon properties of JIG, i.e. 
an artificially electron-doped C2, where the electrons in 
excess are neutralized by a uniform, positive background 
charge; the comparison of the band structure of JIG 
with that of real GICs has shown that two important 



independent variables govern the filling of the interlayer 
band: the interlayer spacing c and the number of elec- 
trons added to the system. In our study, these two pa- 
rameters are represented by the doping n, defined as the 
number of electrons added to C2, and neutralized by a 
jcllium background of the same magnitude (the number 
of C valence electrons is thus 4 + n/2) and the c/a ratio 
of a C2 unit cell with the same in-plane lattice constant 
as CaCg. To minimize the number of parameters we in 
fact kept the in-plane lattice constant fixed throughout 
our study, so varing c is like varying c/a. 



II. PHONONS 

The top left part of Fig. [3] shows the phonon disper- 
sions we calculated for the doping levels n=0 (dotted, 
pure graphite) and n=2/2> (corresponding to CaCg, full 
line). We used the interlayer spacing c and planar lat- 
tice constant a of CaCg in both cases, to compare with 
available experimental and theoretical data. With two 
atoms per cell there are 6 phonon branches, and since 
a phonon displacement is a vector like an electronic Cp 
function, the 6 phonon branches are topologically equiv- 
alent with the electronic C p bands, that is, with the tt 
and tt* bands and with the non-s part of the a and a* 
bands in Fig.[T] In particular, the two phonon branches 
which are the softest along KH, are the out-of-plane vi- 
brations with dispersions similar to those of the tt bands 
in Fig.[T] the four other branches are the in-plane vibra- 
tions. In-plane and out-of-plane modes hardly mix. For 
that reason, the n=0 and n=2/3 phonon dispersions are 
qualitatively similar. 

The only qualitative difference occurs along the T-A- 
L path, where the out-of-plane acoustic branch of JIG 
has imaginary frequencies. Such an instability is an ar- 
tifact of JIG: at n = this branch is just very soft and 
sensitive to computational details, but not unstable^ in 
actual GICS the intercalants contribute to the inter-layer 
binding and stabilize the strucure. Fortunately this ar- 
tificial JIG instability occurs where the e-ph interaction 
is unimportant (no green dots where the frequencies are 
negative in the upper left panel of Fig. [3]), and does not 
affect our subsequent analysis. 

In real CaCg and other GICs 34 in- and out-of-plane C 
modes are also separated in energy, and occur at similar- 
energy scales, and they do not mix with the intercalant 
modes, which are at a lower energy scale (ui ~ 10 meV 
in CaCei 4 -^). 



III. ELECTRON-PHONON COUPLING 

A. Jellium Intercalated Graphite (JIG) 

In a metal, e.g. for doped graphite, the v<\ phonon can 
couple between electrons nk and m (k + q) at the Fermi 
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surface, and thereby attain a relative line-width: 

. nk, rn (k+q)l ^(e«k)<5(emk+q)- (1) 



200 



7^q 



■"'1 

Here, is the average over the Brillouin zone and 

ff^nk,m(k+ q ) = (nk \5V (r)| m (k + q)) /8Q UCL (2) 

is the matrix element of the perturbation, SV (r) , of the 
self-consistent electronic potential due to the displace- 
ment, 
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-SQu 



in the a-direction of the jth atom (and continued as a 
wave with wavevector q), e a j^q being the eigenvector of 
the vc\ phonon. 

The squared e-ph matrix element is usually a smoother 
function of k and k + q than <5(e„k)<5(£mk+q)- In this 
case, it is useful to discuss the so-called nesting functions, 



Xr, 



t(q) = ^28{e nk )8(e mk+q ) 



(3) 
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dk 



e„k=e m (k+q)=0 | v "k X V m(k+q) | BZ 

= - lim n — S x° m (q,^)- 

uj^O 7TLO 

which may be expressed as line integrals (in one 
Brillouin-zone, of volume £Ibz) along the cut of the n 
and the m sheets of the Fermi surface, with the latter be- 
ing displaced by — q. Here, v„k = Ssn^/Sk is the Fermi 
velocity. This is the second line of Eq. [31 Its last line 
relates the nesting functions to the bare susceptibilities. 
Note that the sum over all q simply yields the product of 
the n and m-sheet state densities per spin at the Fermi 
level (e_F=0): 

X»m(q) = 8(e nk )8(e m v) = N n (0) N m (0) , (4) 

q kk' 

because N n (0) = X)k^( £ «k)- Finally, the probability for 
an electron to scatter from and to any point of the Fermi 
surface with total density of states N (0) = ^2 n N n (0), 
due to the interaction with a phonon of frequency w, is 
given by the Eliashberg spectral function 



a 2 F{co) 



(5) 
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^ 8(Snk)8(s m k+q) ^ \9v, nk, m(k+q) \ 2 8(u ~ Uvq)- 



nmlt ^q 

The total e-ph coupling constant is then: 



A = 2 



a 



d0J = A„ q , 



(6) 
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FIG. 3: (Color Online) Upper left: phonon dispersions 
(black solid lines) for JIG with the same lattice constants 
(c/a = 1.81) and filling n —2/3 as in CaC6- Thin dotted 
lines show the corresponding dispersions for n —0 (undoped 
graphite). Imaginary frequencies are shown as negative (see 
text). The area of a colored dot is proportional to the mode 
coupling constant A„ q . Color code: red it* intraband; blue IL 
intraband (negligible); green IL-^k* interband. Upper right: 
JIG e-ph spectral function a 2 F(u>) and X(uj) defined as in 
Eq.|6] but integrated to ui rather than to oo (dashed). Lower 
left: JIG nesting functions Xnm(q). Lower right: Fermi sur- 
face of JIG in the first Brillouin zone (IL-sphere and tt*- 
cylinder) . 



and the last equation uses Eq. [T] for the phonon line- 
widths to express A as a sum over mode coupling con- 
stants 



A 



1 



7^q 



(7) 



l/q 



These quantities are shown for JIG in Fig.(3J in the up- 
per half we focus on the vibrational properties and the 
e-ph coupling (phonon dispersions and partial Eliashberg 
functions) and in the lower half on the electronic proper- 
ties (nesting functions and Fermi surface). 

We see that the Eliashberg function a 2 F(u>) merely has 
two pronunced peaks, respectively centered at w^40 meV 
and w~160meV. The first involves the carbon out-of- 
plane, zone-boundary modes, and the second involves the 
in-plane bond-stretching, zone-center modes. The same 
two peaks (compare our Fig. [3] with Fig. 3 of Ref. [Til ) 
are found in the a 2 F(u>) of CaCg, which, of course, also 
has an additional low-frequency peak involving in-plane 
vibrations of Ca. 

Integrating a 2 F(u)/u} over positive frequencies^ 
yields: A=l.l, thus even larger than in CaCg, where 
A=0.84, despite the absence of the calcium modes. We 
have thus shown that a large e-ph interaction survives 
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in JIG even after removing the intercalants, whereas all 
the recent studies on superconductivity in GICs have as- 
sumed that either the mechanism for superconductivity 
is not conventional e-ph coupling^ or that it is intrinsi- 
cally related to the presence of the intercalanti 13 ' 14 The 
discussion of possible reasons for the quantitative dis- 
crepancy between the A's of JIG and CaCg is postponed 
to the end of the section. 

The Fermi surface of JIG (bottom right panel of Fig. [3j 
has two sheets, both of electron character: the sum of 
their volumes is = VIbz/S- The it* sheet has 

negligible fcj-dispersion and is open, with the shape of 
honeycomb enclosing the vertical KM-HL BZ boundaries. 
Its hollow inside forms a TA-centered cylinder with av- 
erage radius fc^* as shown in red at the lower right of 
Fig. [3j The IL sheet, shown in blue, is three dimensional, 
T-centered, and closed; when n and c/a equal those of 
CaCe, as in Fig. [31 it closely approximates a sphere with 
radius fciL- 

The e-ph interaction can take place in three chan- 
nels: IL intraband, 7r* intraband, and IL-7T* interband. 
The A„ q 's associated with each phonon mode (Eq.[7]) are 
shown in the upper left of Fig.[3J by full dots with area 
proportional to A^q. The contributions from the IL and 
7r* intraband scattering are respectively blue and red, 
consistently with the color code used for these two bands 
in Fig. Q] and for the corresponding Fermi surface sheets 
in the lower right panel of Fig. [3J The interband IL- 
7r* scattering contribution, is, instead, green. The up- 
per right of Fig. [3J shows the Eliashberg function and the 
lower left the nesting functions with the same decompo- 
sition scheme and color code. 

Fig. [3J clearly shows that the zone-boundary out-of- 
plane buckling phonons with w~40 meV give the largest 
contribution to the total A, and that the coupling is 
due to interband IL-7T* coupling. The 7r* intraband cou- 
pling by the optical in-plane bond-stretching phonons 36 
cause the large a 2 F(uo) peak at oj^l60meV. Due to the 
high frequency, however, this peak contributes a mere 
20% to the total A. Finally, with the absence of inter- 
calant modes, our JIG model has no phonons with small 
|q| causing significant coupling (IL intraband and inter- 
band). 

To understand the origin of the large interband e-ph 
matrix element in Eq. [21 we go back to the Wannier-like 
orbital in Fig. [21 As we saw, the IL and tt* bands cross 
essentially without hybridizing (a tiny hybridization oc- 
curs only along the A-L line) because the IL Wannier-like 
orbital is invariant under all operations of the point group 
of the interstitial site, while the sign of the Wannier-like 
function for the tt* band alternates within the extent of 
the IL function. If, however, we apply a buckling distor- 
tion to the graphene sheets, half of the p z orbitals are 
driven closer to the interstitial site where the IL orbital 
of Fig. [2] sits, and the other half are driven away, with 
the net result of a non-zero e-ph matrix element between 
the tt* and IL bands (see the Appendix A, where this is 
demonstrated using a TB model). 



In fact, we only need the e-ph matrix elements for 
Bloch functions at the Fermi level, and not for the en- 
tire IL band. Since the n=2/3 doped electrons not only 
go into the IL band, but also, and even more, into the 
tt* band, the former is less than 1/3 full. To calculate 
the IL Bloch functions at ej? we could therefore also use 
the Wannier-like function for the lowest IL subband, ob- 
tained by folding the IL band into e.g. the three times 
smaller BZ of CaCg. This function is shown in the up- 
per part of Fig. 2] using the same contour as in Fig. [21 
It is seen also to have the full point-group symmetry, 
and since the basis set spanning the lowest IL subband 
consists of the replicas of this function placed at all Ca 
sites, that is, at every third interstitial site, the square of 
this function is three times more diluted (extended) than 
the square of the Wannier-like function for the entire IL 
band. Whereas the Wannier function for the lowest IL 
subband has a\ g symmetry, the ones for the two upper IL 
subbands have e g symmetry, that is, their signs alternate. 
So one could also say that, for these three subbands, the 
Wannier function can be obtained by forming the three 
linear combinations (of axial symmetry a\ g and e g ) of the 
"small" Wannier function of the entire IL band (Fig. [2]) . 
For the lowest band this yields the "large" Wannier func- 
tion shown in Fig, 31 whose axial symmetry is a\ g . The 
point is, that for the relatively low fillings of the IL band 
which are of our interest, this is the Wannier function 
required to form the proper Bloch IL functions at the 
Fermi level; and this function is uniform on the scale 
of the C-C bond length. For that reason there is a large, 
"dormant" IL-7T* matrix element for buckling modes. 




FIG. 4: [Color Online) Top and side view of the Wannier- 
like function for the lowest of the three IL subbands, the ai g 
subband, in JIG (top) and in CaC6 (bottom). The basis set 
spanning the a\ g subband consists of the replicas of this func- 
tion placed at all Ca sites, i.e., at every third interstitial site, 
the blue ones. The isocontour and the color scheme is the 
same as in Fig. [2] 
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FIG. 5: (Color Online) Fermi surface of JIG and of of CaCe, 
both shown in the reciprocal space of rhombohedral CaC6 
(top view, compare to the 3D view of Fig. |3j lower right 
panel). Once folded into the three times smaller BZ, the 
JIG Fermi surface (left), which in the hexagonal BZ (grey 
dotted lines) consists of an IL sphere (blue) and a n* honey- 
comb (red), reveals remarkable similarities with that of CaCe 
(right). 




B. JIG vs. CaC 6 



How well do our results for JIG reproduce what is ob- 
served in real GICs? As an example, we consider CaCg. 
Here, the IL band is folded three times into the rhombo- 
hedral BZ, and split; however, the occupied band struc- 
ture of CaC6 still bares a very close resemblance to that 
of JIG4 ' 13 i 14 ' 17 This is clearly visible in Fig. [71 where we 
show an in-plane cut of the Fermi surface of the two sys- 
tems in the (three times smaller) rhombohedral Brillouin 
zone of CaC6^ The folding of the hexagonal BZ of prim- 
itive graphite into the rhombohedral BZ of CaC6, shown 
in the left panel, may obscure the simple Fermi surface 
topology appearing in the lower right panel of Fig. [31 but 
it clarifies the origin of the complicated Fermi surface of 
CaC6, shown in the right panel of Fig. [7J In both cases 
we clearly recognize a quasi-spherical IL Fermi surface 
and a more complicated structure, which, however, only 
amounts to folding the outer 7r* honeycombs; in the right 
panel only small hybridization gaps, in the T — M and 
equivalent directions, make CaC@ slightly different from 
the purely folded graphite shown at the left. 

This comparison suggests that the introduction of Ca 
atoms between the C2 layers preserves the physical pic- 
ture discussed above. This is also confirmed by the 
Wannier-like function for the lowest IL subband, shown 
in the bottom panels of Fig. [U and is seen to be very 
similar to the one for the lowest subband of JIG, shown 
in the top panels and discussed before. The main differ- 
ence is the presence of some Ca3d3 Z 2_ 1 character in the 
IL Wannier-like function for CaCg, but this does not dis- 
rupt the large dormant IL-7T* matrix element for buckling 
modes. 



C. Why JIG overestimates the CaCe e-ph 

A large e-ph interaction associated to C out-of-plane 
modes has been calculated in CaC6 and several other 
GICs j 14 i 15 ' 34 our JIG model gives a convincing qualita- 
tive picture of this electron-phonon interaction in CaCe 
and, as shown in the next section, also correctly explains 
the observed trend of A observed in different GICs. How- 
ever, as we already pointed out, there is a large quantita- 
tive discrepancy in the values of A, even though the elec- 
tronic parameters (DOS, Fermi velocities) are compara- 
ble. As a matter of fact, in JIG we obtain an overall larger 
A than in the corresponding real GICs, even though, in 
JIG, the phonon modes associated to the in-plane vibra- 
tions of the intercalant, which give a sizable contribu- 
tion to A in real intercalated graphites CaCe, SrCa6, and 
BaCg, are missing (these calcium modes cause a strong 
e-ph interaction, because displacing the IL Wannier-like 
function off center in the xy-plane also awakes the dor- 
mant IL-7T* and IL-IL couplings). 

The bottomline is that, in JIG, the A associated to the 
IL-7T* coupling via out-of-plane phonons is almost three 
times as large as in real GICs. Why? Most of this boost 
comes from the difference in phonon frequency: in real 
compounds, the presence of an intercalant ion partially 
hinders the vertical motion of C atoms, thus hardening 
the out-of-plane frequencies. 

For example, the peak associated to C out-of-plane vi- 
brations is at w~55 meV in CaCe, but at w~40 meV in 
the corresponding JIG. Assuming that the partial A asso- 
ciated to this mode can be approximated by the Hopfield 
formula A = N(Q)g 2 /Mcu 2 , where g is the e-ph matrix ele- 
ment, this would give an enhancement factor of (55/40) 2 , 
i.e., almost two, going from real CaCe to JIG, if N(0) and 
g were the same in the two systems. But they are not, and 
the quantitative differences in the band structure which 
affect them (such as a weak hybridization between IL and 
7r* bands, visible in CaC6 but absent in JIG, see Fig. [Jj 
are probably enough to explain the total enhancement 
factor of almost 3 found between JIG and actual CaCe, 
as far as the e-ph coupling of carbon-buckling modes is 
concerned. 



IV. EFFECTS OF ELECTRON DOPING, n, AND 
INTERLAYER SPACING, c. 

Csanyi et al4 have proposed that the electronic prop- 
erties of donor-intercalated GICs be mainly determined 
by the inter layer spacing and valence of the intercalant; 
recent first-principles calculations 34 and experiments 4 ^ 
suggest a possible correlation between T c and the c axis 
lattice constant; these hints are both confirmed and clar- 
ified in their origin by the results presented here. Having 
so far studied in detail the electron-phonon interaction 
of JIG in comparison with CaCe, we will now turn to 
discuss how the two independent variables of our model, 
c/a and n, affect the e-ph properties of JIG, and compare 



FIG. 6: {Color Online) Position of bottom the IL band in 
JIG as a function of interlayer C spacing (c/a) for different 
electron dopings n; see also Fig. 4 of Ref. @. 



the JIG trends with the available experimental data. In 
the whole range of c/a and n considered here, there is a 
single interlayer band at the Fermi level (see Section I ) . 

Electrons: In Fig. [6] we plot the energy position of the 
bottom of the IL band as a function of c/a for different 
n. As mentioned in Section I, the addition of a uniform, 
positive background of charge to graphite shifts the IL 
band down in energy w.r.t. the rest of the band struc- 
ture, leaving its shape unaffected. The effect is obviously 
stronger for larger n. The IL band is shifted down in 
energy also when the c/a ratio is increased. Changing 
the separation between the graphite planes also sensibly 
affects the IL bandwidth along k z , so that the associated 
FS, which in CaCg is nearly a sphere (Fig. [3]), becomes 
an oblate ellipsoid for smaller c's. In the opposite limit 
(large c's), the opposite happens. The same trend is ob- 
served in real GICs (see Ref. |34{ ). 

Phonons: In the upper left panel of Fig. [3] we show the 
phonon dispersions with electron doping rt=2/3 (black 
solid line) and without doping (n=0, thin dotted line). 
The JIG phonon frequencies are lower than those of un- 
doped graphite in the entire BZ, regardless of the pres- 
ence or absence of e-ph coupling. Such a large, disper- 
sionless energy shift, common to all phonons, is not due 
to e-ph coupling, but rather to a softening of all bonds, 
related to the increased filling of the tt* antibonding 
states^ Fig. [3J shows that on top of this global phonon 
softening, there is a further q-dependent softening caused 
by e-ph interaction. 

The effect of doping and the c/a ratio on the frequency 
of the zone-center (F) optical buckling modes is illus- 
trated in the left panel of Fig. [JJ For fixed c, the mode 
softens for larger n, because of the filling of the antibond- 
ing 7r* band. For fixed n, the frequency softens at lower 
c/a, where the charge transfer to 7r* states increases, due 
to the emptying of the IL band (Fig. [5]). In the right panel 
of the same figure we plot the frequency of the doubly- 
degenerate buckling mode at the (K) point; here, on top 
of the softening due to charge transfer, an additional soft- 
ening due to e-ph interaction takes place. 

The same softening of the out-of-plane phonons with 



1.4 1.6 1.8 2 2.2 1.4 1.6 1.8 2 2.2 

c/a c/a 

FIG. 7: (Color Online) Frequency of the optical buckling 
mode at F (left) and K (right) as a function of the c/a ratio 
for JIG, at different dopings n; the in-plane lattice constant 
is that of CaCs, for which c/a=1.81 and n=2/3. Symbols 
represent calculated data, lines are meant to guide the eye. 
Imaginary frequencies are shown as negative. 



respect to pure graphite is observed also in real GICs, 
where the effect is partially compensated by the presence 
of the intercalant, which hinders the vertical motion of 
the C atoms. However, in some cases the effect depicted 
in Fig. [7] is so strong as to lead to structural instabili- 
ties also in actual compounds. An extreme example in 
this sense is LiC2, for which, using the lattice parameters 
of Ref. 1 381 we found imaginary frequencies in the entire 
BZ for both the acoustical out-of-plane and the optical 
buckling branch. Fig. [7] further suggests that many com- 
pounds with chemical formula AC2 may be dynamically 
unstable for this mechanism. 

Electron-phonon coupling: We have just seen that the 
largest source of e-ph coupling in JIG is the interband IL- 
7r* channel, and that c/a and n critically affect the rele- 
vant electronic bands and phonon modes. To understand 
how this can affect the e-ph properties, we observe that, 
in general, A will benefit from three factors: the filling 
of the IL band (Fig. [5]), the softening of the out-of-plane 
phonon frequencies (Fig.[7J, and an increase of the IL-7T* 
electron-phonon matrix element (<?), which grows as the 
graphene layers are squeezed together (not shown) . This 
suggests that a promising route to increase T c is either 
to increase the electron concentration - but this route is 
severely limited by lattice instabilities - or to reduce the 
c/a ratio, by applying pressure or intercalating smaller 
atoms. Notice however that even in the second route A 
cannot be increased beyond a threshold value, because 
as c/a is reduced the IL band is gradually emptied (see 
Fig. [S]). Experiments have indeed shown that pressure 
can increase T c in these compounds; but this route is 
unfortunately limited by the appearance of instabilities 
in the vibrational spectrum i 15 ' 34 i 39 

To gain a consistent picture of how the interplay of 
the three factors affects the e-ph interaction in JIG, 
we study the Eliashberg function and the frequency- 
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dependent X(u>) as a function of c, with a and n fixed 
at their values in CaCg. In the c-range considered, both 
the position and the shape of the high-frequency peak, 
which involves the tt*-tt* coupled bond-stretching modes, 
is nearly independent of c. For the low-frequency peak, 
involving the IL-7T* coupled buckling modes, the situation 
is entirely different. For the three higher values c/a= 1.81, 
2.00, 2.11 which match those of real GICs (CaC 6 , SrC 6 
and BaCg)^ 1 for which first-principles calculations are 
available^ the contribution to A from the low-frequency 
peak decreases with increasing c: as the graphene layers 
are pulled apart, the overlap between the IL and tt* wave- 
functions decreases, the e-ph matrix-element of Eq.[2]gets 
smaller, the buckling modes get harder. Below c=1.35a, 
on the other hand, the low-frequency peak is just ab- 
sent, because, at n=2/3, the IL band is empty, and gives 
no contribution to A; so at first, as c increases beyond 
1.35a, the e-ph interaction must increase, because, with 
the growing size of the IL sheet, the density of states en- 
tering Eq.[Hgrows. As a consequence A, and thus T c , will 
(for n=2/3) have a maximum for some 1.35 <c/a< 1.81, 
which we calculated around 1.6 and is also shown in 
Fig.[! 
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TABLE I: Comparison of JIG with re=2/3 (this work) with 
real alkali- metal intercalated graphites with rrf/3 (Ref. |34j). 
Eil is the energy of the bottom of the IL band (eV). N (0) 
is the total DOS in states/(eV-C spin). The corresponding 
a 2 F{uj) are in Fig. U 

The data marked with t are taken from Ref. 0; their e- 
ph coupling A also includes intercalant-related modes, which 
contribute almost half the total coupling. 



V. CONCLUSION 



FIG. 8: {Color Online) JIG Eliashberg function a 2 F{uj) (solid 
lines) and frequency-dependent (Eq. |6j but only up to 

ui: dashed lines) for fixed n=2/3 and different c and with a 
fixed at its value in CaC6- For c/a<1.35 the IL band is empty 
and the low-frequency peak associated with the out-of-plane 
C vibrations is absent. As c/a grows from its value (1.81) in 
CaC6 to the one in BaC6 (2.11), this peak shifts to higher 
frequency and its weight dramatically drops, while the high- 
frequency peak of the in-plane modes, remains unchanged. 
Reducing c/a to 1.6, instead, causes an increase in A. 




c/a=1 .35 
c/a=1 .60 
c/a=1 .81 
c/a=2.00 
c/a=2.11 



80 120 

(o(meV) 



200 



In this paper we have studied the e-ph properties of 
electron-doped graphite, using a jellium model in which 
the excess electrons are neutralized by a uniform pos- 
itive background of charge. We have shown that, de- 
spite its simplicity, the proposed model identifies and ex- 
plains an important source of the e-ph interaction in real 
GICs, namely, the existence of an interlayer (IL) three- 
dimensional electronic band, whose filling determines the 
occurrence of a strong, otherwise dormant e-ph interac- 
tion between IL and tt* electrons, for which we provided 
a simple understanding in terms of Wannier-like func- 
tions. Our picture corrects the rough idea according to 
which, since IL states are free-electron like, they must ex- 
perience little or no interaction with the electronic states 
and phonon modes related to the carbon layers. We have 
shown that, on the contrary, IL states can experience a 
strong interaction with the surrounding C lattice, and 
thus give sizable contribution to the e-ph coupling, as 
observed in GICs^ disilicides and metal-sandwich struc- 
tures i 21 ' 22 Our results amount to a natural explanation 
for the empirical correlation between the filling of the 
IL band and the occurrence of superconductivity in the 
GICs, and nicely fit the observed correlation between T c 
and layer spacing in actual GICs. 



Our calculated JIG trend for A matches that of the real 
compounds (see Fig. 7 of Ref. [13). Details are given in 
Table []] which also shows that although this trend is the 
right one, the total A in JIG is consistently larger than 
that of the corresponding GIC. For a discussion of the 
possible reasons for this discrepancy, see the end of the 
previous section. 



Moreover, since tt* and IL electronic bands as well 
as carbon out-of-plane phonon modes are found in sev- 
eral other materials, from nanotubes^ to layered com- 
pounds, 20 including the recently proposed metal sand- 
wich structures ; 21 ' 22 our findings may represent an inspi- 
ration for further research on the electron-phonon mech- 
anism identified here. 
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Appendix A: Tight-binding Model 

Undistorted JIG 

The NMTO method^e a n ows us to define from first 
principles three most localized "atomic" orbitals which 
span the Hilbert space of the electronic states correspond- 
ing to three bands which, in the main paper, were labeled 
as 7r, 7r*, IL, and associated to the e-ph interaction in 
GICs. What do they look like and where do they sit? To 
answer this question we need the geometry of our JIG 
model. "Primitive graphite" C2 has aaa stacking; direct 
and reciprocal lattices are defined by the basis vectors 
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the positions of the two C atoms within the unit cell 
are 



±c ^ ± l/T i _T 2 -T i 
3 V 2 2 



-(. 



±a 



V2 v / 3 



,0,0 



and the position of the interstitial site S, where the in- 
tercalant may sit in real GICs, is: 



s4t, = (o,o,£ 



Out of our three localized orbitals, two sit on the two 
carbon atoms of the C 2 unit cell, are of p z symmetry, and 
are naturally identified with the carbon p z orbitals; the 
third sits on the interstitial site and, although it looks like 
an s orbital (see Fig. [2]) , has only axial symmetry (around 
z); for simplicity of notation, however, we assign it an s 
subscript, keeping in mind that this is not a purely s 
orbital: 




L H K 



FIG. 9: (Color Online) Band structures of undistorted (left) 
and phonon-distorted (right) JIG (see text). Here the energy 
zero is at the Fermi energy. Upper panels: comparison of 
the full band structure based on plane waves and pseudopo- 
tentials (black solid) with the three (tt, n* and IL) minimal 
downfolded bands (red dashed) obtained from a (3x3) NMTO. 
As green dotted lines we also show the corresponding 6 th -n.n. 
tight-binding bands. Lower panels: we reproduce the three 
relevant NMTO bands of the upper panels, and compare them 
with l si -n.n. first-principle tight-binding bands (green solid). 
These bands look obviously poorer than those with 6 neigh- 
bors, but they still correctly reproduce the band splittings 
due to the phonon displacement (see text). 



where the sum extends to all lattice vectors T, we re- 
cover, by construction, the three electronic bands (ir, tt* 
and IL) from which these localized states were derived, 
and obtain the corresponding Bloch states ( | tt; k), \ir*; k) 
and |IL; k)) as a linear combination of these three Bloch 
sums. The three relevant NMTO bands are shown as 
red dashed lines in the upper panels of Fig. [5] for our JIG 
without (left) and with (right) a r-point buckling phonon 
displacement of 0.17 a.u.; the full JIG band structure 
from pseudopotentials and plane waves is shown as black 
solid lines, for comparison. With our 3x3 basis, an essen- 
tialy perfect "first-principles" tight-binding representa- 
tion is obtained by truncating the electronic hamiltonian 
matrix beyond the sixth shell of nearest neighbors (n.n.), 
and shown by the green dotted lines in the upper panels 
of FigO The corresponding formulae and parameters are 
available upon request^ 

However, to understand the basic features of the band 



(r|±C) = ^.(r=FC) 



(r|S)=V s (r-S). 



If we diagonalizc the electronic hamiltonian h in the 3x3 
basis of the localized states | + C), | — C), and |S), or, 
more precisely, of the three corresponding Bloch sums 



i±c,k) = J2 

T 

|S,k) = £ 



Jk-T 



T±C) 



e lk - T |T + S) 



(8) 
(9) 



TABLE II: 3x3 nearest-neighbor first-principles TB pa- 
rameters for undistorted JIG, from NMTO (see text and 
Table |Tn}, in eV. Unlike Fig. d here the energy zero is at 
the crossing point of the tt and tt* bands, so Ef = 2.0 eV. A 
carbon-buckling phonon of amplitude u — 0.02 a (where a is 
the planar lattice constant of C2) changes t sp by as much as 
0.1122 eV. To linear order this means 5t sp /5u = 2.241 eV/A. 



3.5115 
-0.8811 



0.8155 
2.5626 
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TABLE III: The 3x3 electronic hamiltonian for JIG when only nearest neighbor integrals (see Table [TT] and text) are retained. 

|S;k> |C;k) |-C;k) 

(S; k| (s s -e p ) + 29 s [cos(2A) + cos(2^) + cos(2i>)] -2ie~^t ap sin(£) (l + e - 2l "+ e' 2 ^) -2ie~ l h ap sin(£) (l + e" 21 ^ e" 2lA ) 

(C;k| 2ie^t sp sm(0{l+e 2lfl + e 2lv ) - t pp (l + e 2lA + e~ 2iv ) 

(-C;k| 2ie l h ap sm(0(l + e 21 " + e 2lA ) -t pp (l + + e 2iv ) 



structure and e-ph interaction, nearest-neighbor matrix 
elements may be enough. In this case only a few integrals 
survive: if a and /3 equal s for the |S) orbital which sit 
on the interstitial site and p for the two carbon | ± C) 
orbitals, we have onsite hopping integrals indicated by e a 
and £/3, hopping integrals between orbitals of types a, (3 
sitting on n.n. atoms indicated by t a p (and all taken 
positive, t = \t\), and hopping integrals between orbitals 
of type a sitting on the same atom of n.n. unit cells 
indicated by a . Their numerical values are in Table HT1 
Then the hamiltonian matrix is given by Tabic IIII1 
where 



with eigenvectors 



T a 

A = k • — = - (k y + V3k x ) , 



k-^-^k 
K 2 ~ 2 y ' 



v = k- 



-Ti 



[i - A = - (ky-^/Zkx) 



the corresponding n.n. bands are shown as solid green 
lines in the lower panels of Fig. [51 where they are com- 
pared to the original NMTO bands (red dashed) . Notice 
that, in this model, the k z dispersion entirely comes from 
the s—p hopping between the s orbital sitting on the in- 
terstitial site and the carbon p z orbitals on the graphite 
planes immediately above and below it: carbon p—p ma- 
trix elements between adjacent graphite planes are be- 
yond n.n. In the k z =0 plane (£=0) the matrix element 
between s and p z orbitals is identically zero, because the 
phases of the p z orbitals on two adjacent graphene sheets 
cancel. Therefore, for k z =0, the IL band is 



e 1L {k)=e s + 26 s 



k y a + ^/3k x a 



cos k y a 



k v a — y/3k x a 



, (io) 



and its eigenvector |IL; k) = |S; k) just coincides with the 
Bloch sum of Eq. [9j the 7r, 7r* bands entirely derive from 
the p z sub-block, whose diagonalization yields: 



e 7r* ,7r (k) — £71- =p (k) — Ep ± tpp 1 1 + e + 



(11) 



— £p i tpp \ 



kqiCL 



V~3k x 



2 cos k y a 



|7T*;k) = |7T_;k) = i={|-C;k) e ^( k )-|C;k) e -^ k )} 
k ;k) = K + ;k) = -L{hC;k)e^ k ) + |C;k) e -^ k )}, 



where 2ip (k) is the phase of (l + e 2iA + e~ 2lu ) . The state 
|-/T + ) has lower eigenvalue and is of bonding character, be- 
ing a bonding combination of the carbon p z orbitals at 
k = 0; the other one, 1 7r } , is the antibonding combina- 
tion. In Fig. [TU] we show a real-space cartoon of these 
three complex Bloch sums, only in the xy plane; in the 
vertical direction (not shown) both Bloch sums |7r±;k) 
are odd with respect to the plane, since both derive from 
p z states; orbitals on different planes differ by a phase 
factor e 2 ^. 

For k z 7^ the matrix element between s and p z 
states is not zero, so, unlike the k z — case, the |S;k) 
and |7r±;k) Bloch sums may hybridize. A priori this 
means that the matrix elements of the electronic hamil- 
tonian (ir±; k|/i|S; k) may be nonzero; however it turns 
out that one of the two, the antibonding Bloch sum, gives 
(7t_; k|ft,|S; k) = almost everywhere. The only excep- 
tion is a tiny k-space region between the A and L points; 



e itf(k)+2i p _|- e -!f (k)+2iji 

■UpQk) ± e ~" plk) 



FIG. 10: (Color Online) Schematic real-space picture of the 
three Bloch sums |S; k) and \n±; k) in a single graphitic plane 
(see text). Red dots mark the position of carbon atoms; the 
blue dotted circle represents the interstitial localized orbital 
|S) sitting in the central unit cell; similar wavefunctions are 
at the center of all other unit cells, multiplied by an appro- 
priate phase factor, but are not shown here. In the vertical 
(out-of-plane) direction, wavefunctions are multiplied by an 
additional phase factor e 2 '^. 
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FIG. 11: (Color Online) Schematic real-space picture, in the 
xy plane, of the Bloch sums |S;k), |7r+;k) (left), and |7r_;k) 
(right), for k x = k y = -> <p(k) = (T - A line). 

the corresponding, small gap between the tt* and the 
IL bands is visible in Fig. [5] but not in Fig. [TJ because 
of the thickness of the "fat bands". The reason why 
(7r_; k|/i|S; k) =0 may be immediately understood along 
the T — A line (k x = k y = 0), where v?(k) = (see Fig. ITU 
where we display the three corresponding Bloch sums): 
each S site has 12 carbon nearest neighbors, 6 above and 
6 below, which respectively carry, in the vertical direc- 
tion, phase factors of 1 and e 2 ^, so for the bonding (+) 
and antibonding (— ) Bloch sums we obtain: 

(7r+;k|ft|S;k> = -4= (Qt sp + 6t sp e 2 ^) 
v2 

(tt_; k|A|S; k) = --= [3t sp - 3t sp + (3t sp - 3t sp )e 2li ] = 
v2 

However, even for any other k point in the Brillouin zone 
(BZ) a similar, just more lengthy procedure, based on the 
phases of Fig. [TU1 gives zero for the antibonding Bloch 
sum; except, as mentioned, in a small region between 
the A and L points, where the definition of bonding and 
antibonding looses its meaning. In other words, the elec- 
tronic states corresponding to the tt* band may legiti- 
mately be identified with the |7r_;k) Bloch sum of p z 
orbitals over almost the entire BZ, with no or very little 
admixture of the Bloch sum of the interstitial orbitals 
|S;k). 

Distorted JIG 

We have seen that, in the undistorted lattice, the ma- 
trix element of the electronic hamiltonian between the 
Bloch sum of interstitial orbitals |S; k) and the antibond- 
ing Bloch sum |7r_;k) is almost always zero, so that the 
latter can be practically identified as the Bloch state of 
the 7T*, while the Bloch state of the IL band will have 
a prevailing |S;k) character plus, at certain k's, some 
admixture of the bonding |7r + ;k) Bloch sum. Upon out- 
of-plane distortions of the graphene sheets, our e-ph ab- 
initio calculations reveal, instead (see main paper), a 
large coupling between the the tt* band and the IL band. 




FIG. 12: (Color Online) Schematic picture, in the xy plane, 
of the Bloch sums |S; k), |7r+; k) (left), and \tt-\ k) (right), at 
the r point of graphene with an optical buckling distortion of 
q=0. 

How do we understand that? By symmetry: in graphene 
the relevant out-of-plane phonons have the same sym- 
metry as the tt* electronic states. The relevance of this 
symmetry argument appears if we estimate the e-ph cou- 
pling between the electronic states i and j due to the 
phonon mode (f q) with the change of the corresponding 
electronic matrix element upon phonon displacement 

(i,k\5h ucl \j,k + q) , 

multiplied by the nesting function <5(£ik)<5(£jk+q) and in- 
tegrated over the Fermi surface. Within our nearest- 
neighbor model this well defined, but not immediate pro- 
cedure, ultimately recovers the high coupling obtained in 
the full first-principles calculation (see main paper) . But 
at special points in the BZ our model allows to under- 
stand the underlying symmetry argument by inspection. 
The simplest example is the T optical buckling phonon 
(q = 0). This is the frozen- in distortion considered in the 
right panels of Fig. [51 whose most visible effect on the 
electronic states is the appearance of a hybridization gap 
between the tt* and IL bands at ~ 2 eV above the Fermi 
energy, in the middle of the T — K and T — M lines. If, 
upon a q = phonon distortion, the electronic hamilto- 
nian is written as h u = h u= o + 5h u , then for a general 
k^O the change in the electronic matrix element relevant 
to the e-ph coupling becomes (i, k| 8h u \j, k) and is given 
in Table [TV] From it one may already get some intuition 
of the strong coupling between the tt* band and the IL 
band. 

At k = the symmetry argument becomes obvious: 
the induced by the optical buckling phonon (al- 

ternate raising and lowering of carbon atoms) has the 
same symmetry as the |7r*;k = 0) = |7r_;k = 0) electronic 
state (alternate sign on carbon atoms, antibonding com- 
bination of p z orbitals); the product of two anti-bonding 
functions is bonding; as a result, we have a strong over- 
lap with the IL electronic state, which, at k = 0, is the 
simple sum of interstitial |S) orbitals. This is schemati- 
cally shown in Fig 1121 The only difference with respect 
to the purely electronic analysis of the previous cartoons 



12 



|S;k) |C;k) |-C;k) 

<S;k| 2e- I «^ sp cos(£)(l + e- 2v +e- 2 « / ) ~2e~^5t sp cos(£) (l + e" 2v + e" 2jA ) 

(C;k| 2e^<5t sp cos(£)(l + l + e 2l " + e 2 " / ) 

(-C;k| -2e l «^ sp cos(£) (l + e 2lfl + e 2lA ) 



TABLE IV: Change in the JIG nearest-neighbor tight-binding electronic hamiltonian matrix elements under a F-buckling- 
phonon distortion (see text and Tables ILTl and [TTT)) . 



(Figs. [TUl ITT|) is that, besides the electronic phases, we 
must now consider the phonon phases too. Here St sp in- 
dicates the change of the corresponding matrix element 
upon distortion; quantitatively, for a distortion ampli- 
tude 
St. 



0.17 a.u., the t sp matrix element changes by 



1 Sp 



<>sp 

0.101 eV. For those carbon atoms which move 
towards the interstitial site, where |S) sits, the matrix el- 
ement becomes more negative: —t sp — > —(t sp + 5t sp ); for 
those which move away from it, it becomes less negative: 



-t 



sp 



sp 



St 



sp) 



We then have, respectively for the 



bonding and antibonding Bloch sums: 

(7r+;k=0|5/i„|S;k=0) = 

-^[l + (-l) + l + (-l) + l + (-l) 

(7r_;k = 0|(5/i„|S;k = 0> = 



2St 



V2 



sp [1 - (-1) + 1 - (-1) + 1 - (-1)] = ~^=5t sp . 
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In conclusion, with our simple TB we obtain, at first 
glance for the T point, and with some additional work at 
any point in the BZ, a simple physical picture not only 
for the electronic states, as seen in Figs. [10] and [Til but 
also for the e-ph interaction. 



For the self-consistent calculation of the band struc- 
ture, Fermi surfaces, phonon dispersions and e-ph cou- 
pling 24 of jcllium-doped graphite we employed ultra-soft 
carbon pseudopotentialsj 4 ^ with a cutoff or 30 (360) Ryd 
for the wave-functions (density) respectively, and PBE 
exchange-correlation functional^ For the k-space inte- 
gration of the self-consistent calculations we employed a 
16 x 16 x 8 Monkhorst-Pack grid, 45 yielding 150 points in 
the IBZ, with a 0.06 Ryd cold smearing-^ a much denser 
32 x 32 x 16 grid was used for the integration on the Fermi 
surface of the e-ph matrix elements. This led to a conver- 
gence of ~1 meV for the phonon frequencies, except for 
the acoustical out-of-plane branch along the T — A line, 
which requires a much higher cutoff to converge. For con- 
sistency check, we also calculated the T— point frequen- 
cies of graphite in the experimental structure, obtaining 
frequencies within ±0. 1 meV of Ref [33j. Dynamical matri- 
ces and phonon linewidths were calculated on a 16 x 16 x 8 
grid in q-space and the phonon dispersions were then ob- 
tained by Fourier interpolation. Such a dense mesh in 
q-space was necessary to resolve the effect of the Fermi 
surface nesting on the e-ph coupling (Fig. [3]). 



Appendix B: Computational Details 

The results presented here were obtained within den- 
sity functional theory^ using two different ab-initio 
methods: for "fat bands", Wannier-like orbitals, and 
tight-binding parameters, we used TB-LMTO-ASA 25 
and NMTOr& and for total energies, phonon frequen- 
cies and e-ph couplings, we employed plane- waves and 
pseudopotentials^ All methods give the same results for 
the electronic bands, the Fermi surface, and several other 
properties, as was carefully cross-checked. 



For the fat bands and Wannier-like orbitals, we used 
the TB-LMTO and NMTO methods^^ l n this case, 
to obtain a good representation of the electronic wave- 
function both in the neighborhood of C atoms and in the 
interstitial regions, we had, in our C2 unit cell, an s,p, d 
basis set on each C atom, but also two "empty spheres" 
(with s,p channels): one, small, in the C plane; another, 
very large, between the planes, on the interstitial site 
called S in appendix A, whose crystal coordinates are 
(0,0, |). This is where the intercalant may sit in real 
GICs: for example, on each of these interstitial sites in 
LiC2; on every third such sites in CaCg. 
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